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I. INTRODUCTION 



The possibility of producing mini black holes in high energy collisions is one of the exciting 
prospects of research at the LHC. It appears as one of the major implications of brane- world 
models, where thepresence of extra dimensions opens up the possibility for a lowered Planck 
scale 

mm a], so that gravity is much stronger in the spacetime which includes the 
extra dimensions and black holes can then be produced in high energy collisions [6 



Two key signatures of such black holes would be the observation of Hawking radiation 
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3,[i3l,Q 



Although at 



and the measurement of quasi- normal mode frequencies 1^, 
the horizon of a black hole the emitted Hawking radiation is a perfect black-body radiation, 
for an observer at infinity the radiation appears modified. This is because in escaping from 
the black hole the radiation needs to cross the non-trivial potential around the black hole 
which implies that part of the radiation will backscatter, resulting in a frequency-dependent 
factor that modifies the otherwise ideal black body radiation. This is known as a grey-body 
radiation and has a spectrum described by a modified black-body emission rate, which in 
D = 4 + n dimensions is of the form 

dE(u) ^ , , CO d"+3fc 



e'^/Tn ^ 1 (27r) 



n+3 



(1) 



where i stands for the angular momentum quantum number, and k for any other quantum 
numbers of the emitted particles. ct£^k(i^) is the grey-body factor modifying the black-body 
radiation at the horizon to the spectrum observed at asymptotic infinity away from the black 
hole. The Hawking temperature Th for a Schwarzschild black-hole is given in terms of the 
horizon radius Vh and the number of extra dimensions n 1^ , 

n + 1 



Anrh 



(2) 



n + l 



(3) 



For such a black hole, Vh is enhanced and given by 

^ 1 /MBH\^/^8r(^) 

y?M, \ M, J y n + 2 

where Mbh is the mass of the black-hole. 

Clearly, an experimental observation of Hawking radiation from such black holes is a very 
exciting prospect as it constitutes direct observational evidence of black-holes. It should 
provide a rich source of information about the dimensionality and structure of spacetime as 



well as constitute a first step toward establishing quantum gravity as a discipline accessible to 
experimental investigation. One way of exploring experimental implications of such theories 
that can be observed at accelerators has been through event generators. Many such programs 
have appeared over thevears fl H Q. the most rece,. a„d compreheusn-e of which 
is perhaps BlackMax [IS]. Based on phenomenologically realistic models free of serious 
problems that plague low-scale gravity, it incorporates all grey-body factors known to date 
along with several effects that make predictions which are as realistic as possible. 

In brane-world models, we need to localize Standard Model (SM) particles on a 3-brane 
yet allow gravity full access to the bulk. Such a localization is often done via a projection 
from higher dimensional spacetime onto a four-dimensional hypersurface identified as a 3- 
brane representing our world. While this method succeeds in describing the geometry, it 
fails to capture the defining feature of the brane, which is its tension. The latter is often 
neglected owing to the difficulties in obtaining exact solutions of Einstein's equations in 
higher dimensions in the presence of a brane of finite tension. Recently, a solution that 
incorporates a 3-brane of finite tension in six-dimensional spacetime was constructed by 
Kaloper and Kiley for static black holes 20||. It was subsequently extended to include 



rotating black holes 



2l| . In this model, at the position of the brane a conical singularity 



forms that acts to off-load tension thereby keeping both bulk and brane locally fiat. Some 



work 



22 



23 



24 



25 



26| has already been done to study the implications of including tension 
in this way, resulting in qualitative as well as quantitative modifications which may affect 
the interpretation of experimental results. 

The present work aims at continuing in that line by studying analytically the effects of 
including tension on grey-body factors. In detail. Section [IT] is a quick review of the main 
effects of including tension in the Schwarzschild case and the attendant modifications of 
the parameters of the black-hole. In Section IIIII we obtain analytical expressions for low- 
frequency grey-body factors for bulk as well as brane-localized emissions, and and then 
compare with exact numerical results. In Section [IV] we calculate analytically grey-body 
factors in the case of large imaginary frequencies, again for emissions both in the bulk and 
on the brane. In the tensionless limit, our expressions are valid for arbitrary number of extra 
dimensions. Finally, Section |V] contains our conclusions. 
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II. SCHWARZSCHILD BLACK HOLE ON A TENSE CODIMENSION-2 3- 
BRANE 



Black-holes formed through high energy collisions are expected to evaporate through 
several phases (e.g. P]). After the initial "balding" phase where the gauge hair and other 
asymmetries are lost, a spin-down phase ensues, with a small portion of mass being lost 
by radiating away the angular momentum. However, it is expected that the black-hole will 
spend most of its life in the Schwarzschild phase emitting spherically symmetric Hawking 
radiation. In what follows we will concentrate only on this phase. 

In brane-world scenarios, a small black hole with radius which is much smaller than the 
size of the extra dimensions R is in fact completely immersed in the full 4 + n dimensional 
spacetime. In asymptotically flat space, the metric that describes a non-rotating black-hole 



id, 



in higher dimensions is a generalization of the 4-dimensional Schwarzschild metric 

d.^ = -Ur)dt' + -l^dr^ + r'dQl^^, /„(r) = 1 - (^) (4) 

where the element dVL2+n defines a (2 + r;,)-sphere. To geometrically localize fields on the 
brane, what has often been done is to project onto the brane by setting all angular coordi- 



nates 6i to 7r/2, except for the usual 4-dimensional angles 6 and ip (see, e.g., |27|, l28|). The 
resulting 4D metric then reads 

ds' = -Ur)dt' + j^dr^ + rMfi^ (5) 

As mentioned earlier, such treatments neglecting the tension A of the brane are mainly due 
to the the difficulty of obtaining higher-dimensional solutions of Einstein's equations in the 
presence of a brane of finite tension. They can be justified on the grounds that in order for 
a black hole to be treated semi-classically, its mass must be much bigger than the Planck 
scale, so that the tension of the brane does not appreciably alter the black hole background 
and at a scale 0{rh) the solution is fiat. While this is in principle true, it turns out that 
for a black hole to be clear of the quantum regime in this context, given a gravity scale 
of the order of 1 TeV, the mass of the black hole needs to be on the order of just a few 



TeV 



However, because brane tension will be generated from contributions to vacuum 
fluctuations due to brane-localized fields, one expects the tension to be of the same order as 
the fundamental scale of the theory. Therefore, in order to be realistic, the brane tension 
ought to be taken into account. 



By introducing tension using the Kaloper-Kiley codimension-2 setup [20|, the tension 
of the brane does not curve either the brane or the bulk, but forms a conical singularity 
instead, creating a deficit angle Aip. By adapting to brane tension, Aip provides a self-tuning 
mechanism which leaves both bulk and brane locally flat. It is given by 

A^ = 27r(l-6)AM,^ (6) 

where < 6 < 1. The resulting geometry is the same as that of Schwarzschild, except that 
a wedge defining the deficit angle is now cut out with its edges identified, or equivalently, 
the coordinate gets rescaled. The resulting metric reads 

dsl = -f{r)de + f{r)-V + r^dnl /(r) = 1 - (^) ^ (7) 

with 

dnl = de^ + sin^ e [d^p^ + sin^ ^{dx^ + sin^ xd^^)] (8) 

One significant implication that can be readily shown is that the radius of the horizon gets 
enhanced, 

hence also the geometrical cross-section and production rate, which is important from an 
experimental point of view. 

Semi-classically, we can study the propagation of particles of spin s in the vicinity of a 
black hole by writing the corresponding field equation in the background of the metric of 
the black hole. The resulting equation can then be separated into a radial equation and a 
set of angular equations. The radial equation is most conveniently dealt with by casting it 
into a Schrodinger-like form, 

+(u;'-V [r (r,)]) ^ (r,) = (10) 
where is the "tortoise coordinate" defined by 

(11) 



/(r) 

and V{r) is an effective potential which depends on the type of particle considered. The 
potential also depends on the eigenvalues of the angular equations, so that in order to find 
the observationally significant eigenvalue u of the radial equation, we need to first solve the 



angular eigenvalue problems. An exact solution for the angular eigenvalues for integer spin 



was first obtained in 
fermions in 30|. 



29|. It was subsequently derived using Jacobi polynomials for Dirac 



We can sum up the effects of including tension in the codimension-2 model for emissions 
in the bulk by the following replacement rules 

r,^^, i,^\, = i,+ (l-l]m, ij^b^P (12) 



^' ' \b 

where ii is the orbital quantum number for the ith angular equation, and m is the magnetic 
quantum number. For brane-localized emissions the same rules apply except that the orbital 
quantum number does not get modified. In subsequent sections these rules will be the main 
ingredient in obtaining explicit analytical expressions for grey-body factors including brane 
tension. 



III. LOW-FREQUENCY GREY-BODY FACTORS 



Calculations of grey-body in the case of low-frequency have been performed both for bulk 



271 . l3ll . |32| and brane-localized emissions 



28 



33, 



34| . In this section we shall extend these 



results to include brane tension. 

In general, one tries to find an approximate solution of the radial equation which is valid 
at infinity and another one valid near the horizon. Then by matching the two solutions in 
the common domain of validity one can extract the coefficients of reflection and transmission 
and hence deduce the absorption coefficient (grey-body factor) as a function of frequency. 
For bulk emissions of integer spin, we can use the Ishibashi-Kodama 'master equation' 



321 . |35|, in which the potential in D = n + 4 dimensions is given by 



V = f{r) 



+ L>-3) {D-2){D-A) (l-j92)(D-2)2r 



2 ^D-Z 



+ 



+ 



47-2 ' 4 

with p determining the type of field (particle) one is interested in 

, scalar and gravitational tensor 

2 , gravitational vector 

P= \ 

2/{D - 2) , EM (gauge) vector 

^ 2{D - 3)/{D - 2) , EM (scalar) 



„D-1 



(13) 



(14) 
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and ^ being the angular momentum quantum number taking on values 



0, 1, . . . , scalar 



1,2, 



2,3, 



EM 

gravitational 



(15) 



Gravitational scalar fields are governed by a more complicated potential, but they can be 
described approximately using the above potential by taking p to be [32| 



Pgrav-scalar ~ 2 + 0.674D-0-5^^5 



(16) 



Specifically, for the case we are interested in, n = 2 we have 

Pgrav-scalar 2.25 (l"^) 

Writing the radial equation in terms of R{r) = r~^~tvl/ (r^) for n = 2, we end up with 



A3 (3 + A3) 



9 1 



R{r) = 



(18) 



where we have replaced i by A3 (eq. (fT2|) ) in order to account for the brane tension [29|]. Near 
the horizon, it is convenient to make the change of variables r — /(r) (note that /(r) 
as r — > r/j). The resulting equation is 



d^ d 



A3 (3 + A3) , 4 , ^ ^ 



whose general solution can be written in terms of hypergeometric functions, 
i?NH(/) = A_r (1 - ffF{a, b, c; /) + A+r-(l - ffFia - c + 1, 6 - c + 1, 2 - c; /) (20) 



where 



a 



a = a + 13 - -p 



a + l3+-p 



. , . A3 + - 

2 3 y V 2 

3 2 3 



"3 2 3 



A3 + ^ j - {^^Thf - 



A3 + 



(21) 



c = l + 2a = l + -icur,, 
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Imposing the boundary condition that at the horizon the solution be purely in-going, we 
obtain the constraint 

A+ = (22) 
Far from the horizon (/ — 1), the asymptotic solution fl20p under the constraint fl22]) behaves 



as 



rfe\3/3 r(l + 2a)r(l-2/3) 

r) r(l + a-/3 + |p)r(l + a-/?-|p) 



^ 'rv.\ 3(1-/3) r(l + 2a)r(2/3-l) 



(23) 



r/ T{a + (3 - p)T{a + (3+ jp)_ 

which follows from standard identities of hypergeometric functions. 

In the far field region, as r — > oo, or correspondingly / 1, the radial wave equation 
may be approximated by 



9"{r) + -g'{r) + 
r 



-^-(^ + A3(A3 + 3))i 



gir) = (24) 



where we defined g{r) = r^^'^R(r). The solution can be written in terms of Bessel functions, 
^FF(r) = ^ = (5+ J,^+3/2(a;r) + 5_n3+3/2(^r)) (25) 
In the low frequency regime, 

UTh < 1 (26) 



the asymptotic solution ( 12511 may be approximated near the horizon by 

, A3+3/2 



(I) .M « ^(A, + 3/2)p^-'-'''^.-A.-3 



Matching the two asymptotic solutions in the intermediate region ( fl23|) and fl27|) . respec- 
tively), we obtain 

B+ _ r(l - 2p)T {a + p-lp)T{a + f3 + lp)T (A3 + 5/2) T (A3 + 3/2) / 2 ^ 



B- vrr(l + a-/3+|p)r(l + a-/5-|p)r(2/3-l) \ujrh 

(28) 

where we used (3 ~ — A3/3 (eq. ([21])), on account of (^E^. 
The reflection coefficient is 

outgoing amplitude 5+ — iB_ 
incoming amplitude + 
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from which we deduce the bulk absorption probabihty 



bulk! 



\n\ 



2i {B* - B) 



(30) 



BB* +i{B* -B) + l ' 
After some algebraic manipulations, in the low frequency regime fl26l) we arrive at the explicit 

expression for the absorption probability, 



A 



r(i 



A3+2p 
3 



)r(i + ^) 



r(i + ^)r(A3 + |l 



27 



(31) 



31 



36|). 



which reduces to the tensionless result in the limit 6-^1 (see, e.g.. 

Evidently, the dependence of the absorption probability on the brane tension is compli- 
cated. To study its behavior, let us first look at some graphical representations of this result. 
Figure [T] depicts the absorption probability as a function of b with uj fixed. As can be seen, 
all types of perturbation exhibit the same qualitative behavior with the probability vanish- 
ing as 6 — s> and increasing monotonically as 6 1. This behavior persists qualitatively 
at high frequency as demonstrated in figure [2] which was obtained from points generated 
by exactly solving for absorption probability numerically. A further demonstration of how 
tension suppresses the absorption probability is provided in figure [3l where probability vs 
frequency in the case of a scalar perturbation is plotted for various values of b in the small 
frequency regime. The graph also has points representing exact numerical calculations. As 
expected, the predictions of the analytical formula diverge from the exact values with in- 
creasing frequency. Finally, figure H] shows the effect of the magnetic quantum number m 
in the case of scalar emission. The qualitative behavior of all curves is the same, however 
the absorption probability varies over many orders of magnitude depending on m. In figs. 
[3] and HI the type of perturbation chosen has a negligible effect on the outcome. 

Moving on to the case of brane localized modes (the photon, massless fermion and scalar), 
we note that because we are projecting onto the brane, angles x and ip "will be fixed, so that 
the dependence of the metric on b will only be through inside the metric function /. Aside 
from that, the results pertaining to the grey-body factors will be identical to those in the 
tensionless case. To find these grey-body factors we follow the same procedure as above, but 
instead using an equation analogous to the Teukolsky equation on the brane. For a static 



background, the radial wave equation is (see, e.g., |28l. l37l|) 



(/r 



2\l-s 



dRs 
dr 



+ 



f 



2isujr 



isuor 



f 



(n + l)(l-/) 



(32) 



-A-(2s-l)(s-l)(n + l)(l-/)}i?, 
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where A 



d 



Sis 



sin 9 de 



smO- 



d9 



1), and s is the spin. The corresponding angular equation reads 

2 



+ 



2ms cot 9 



m 



smi 



+ s- s^ cot^ 6 + A 



sm 



nm 







(33) 



where e'""^ sS^id) = sl"r (^2) are 
s = h, the absorption probabihty is 



mown as the spin-weighted spherical harmonics. For 



28j: 



U 



brane I ^ 
=1/2 I 



An 



and for s 



I A brane I 



9 



V25/3y [r(j + i)]2 

,2.>2 r(t)r(^)r(j- 
r (^) r (2j + 2) 



2) 



(34) 



(35) 



where j is the total angular momentum in both cases. 
For brane- localized scalars the result is 27l |: 



I .brane 1 2 _ 16vr fUjrh\^j+^ 

' " 9 I 2 J 



^s=0 

Upon comparison of fl5S]l and fl5^ . we obtain 

I /) brane I ^ 



r(^)r(i + t) 
r(i+,)r(i + ^) 



(36) 



brane I ^ 



s=l 



1+J 



(37) 



which shows that (in this frequency regime) the gauge vector and scalar behave the same 
way up to a multiplicative factor, which tends to unity for large j. One can see from the 
expressions above that only the first few values of j contribute to the emission spectrum 
making the emission rate for scalars higher than for gauge vectors in this frequency range. 
This can be seen numerically in the absence of tension in the results of [38]. At higher fre- 
quencies the situation reverses, and the rate for scalars becomes smaller. Graphs analogous 
to those for bulk emissions are shown in figures [5] and [6] for scalars. For low frequencies, 
the behavior with varying tension is now reversed and emission is enhanced with increased 
tension. For higher frequencies the trend is not as conclusive. The difference between the 
two cases can be understood in terms of two competing effects: (a) the enhancement of the 
geometrical cross-section with tension implies an increased rate of emission for both bulk 
and brane- localized modes, and (b) the fact that the angular eigenvalues increase with ten- 
sion and appear only in the bulk potential (causing the potential barrier to increase) results 
in suppression of only bulk emissions. 
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From an experimental point of view, it is of great interest to know the relative emissivities 
between brane and bulk modes. The above discussion seems to hint that the inclusion of 
tension favors the EHM conjecture jsol, namely, that black holes emit mainly on the brane. 
To verify this in our case, let us consider the emission spectrum as defined in eq. ([T]) for 
scalars. For that, we need to write the grey-body factor ai^^lu) in terms of the absorption 
probability and scalar multiplicities Ni^m- For either bulk or brane-localized emissions 
it is of the form 

dtdu 27r ^ e'^/^H _ i 

In the bulk, the scalar multiplicities are |Q] 



m 



0, i (4 + 2) (4 + 1) 



N!:t={ ' (39) 

m ^ 0, (4 - m + 2) (4 - m + 1) 

whereas on the brane they they are equal to the standard value 2i+l. Starting with the bulk, 
we numerically solve eq. (ITS!) for the absorption probability with p = 0, and using eq. (!39ll . 
we end up with the spectrum shown in fig. [8l A similar calculation for brane-localized 
scalars is done by projecting the bulk equation on the brane, and the result is shown in 
fig. [HI Notice that in both figures, the emission rates are enhanced with tension. The reason 
for this comes primarily not from the grey-body factor, but rather from the denominator in 
the black-body distribution which depends on the Hawking temperature that is common to 
both types of emission. We can discount this effect by looking at the relative emissivity, in 
which we take the ratio of the emission rate in the bulk to that on the brane. Figure [TOl 
shows this for scalars. 

Similar plots are obtained for gauge vectors. On the brane, we again have the multiplic- 
ities 2i + 1, but in the bulk we have 

f m = 0, Us (Sis + 5) 

m ^ 0, 3 (4 - m + 2) (4 - m + 1) 

In the bulk, we solve eq. fllSp with p = |, and on the brane eq. (I32p is solved with s = 1. 
The results of doing this are shown in figures [TT] through [131 In this case the emissivity is 
only partial, since the scalar degree of freedom is not included. 

We can see that in both cases, in the range shown, the emission rate is greater on the 
brane than in the bulk, in support of the EHM claim. It should be pointed out, however, that 
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even though mini black-holes at the LHC will spend most of their life in the Schwarzschild 
phase, they first form with non-zero angular momentum. This has been shown to modify 
the Hawking radiation due to super-radiance in such a way that emission in the bulk at the 



initial phase dominates over emission on the brane 



40 



41 



42|. 



IV. HIGH-FREQUENCY GREY-BODY FACTORS 



Next we turn our attention to grey-body factors at large imaginary frequencies. We shall 
discuss different values of spin and take into account the extra dimensions as well as the 



brane tension in the codimension-2 model 
Andersson and Howls 



20|. The calculation will follow the method of 



43| which was then used by Cho [4^. It is a combination of the 



45l | and the standard complex-coordinate WKB 



monodromy argument of Motl and Neitzke 
method. 

The first step is to rewrite the master equation ( ITOi) in a form suitable for the application 
of the WKB method. To this end, we introduce the redefinition of the wavefunction. 



which reduces ( TTUl) to 



The resulting WKB solution is 



0"(r) + g2(r)0(r) =0 



1 



exp 



±i I Q (r') dr' 



(41) 



(42) 



(43) 



with t being a reference point, and 



Integer spin particle are described 



Following Motl and Neitzke 



-V{r) + \{f{r)f -\f{r)f'\r) 



Dy the potential (fT3l) . 



(44) 



45], we go around a contour in the complex r-plane and 



impose preservation of monodromy. Because of the exponential nature of the WKB solutions, 
dominance will be exchanged between the two exponentials as we move around the complex 
plane. This implies that some terms that are exponentially small and can be overlooked 
in one region of the complex plane can grow exponentially in a different region. To take 
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this into account, one needs to resort to the Stokes phenomenon which keeps track of these 
changes. 

In such an analysis, we need to pay attention to the zeroes and poles of Q{r). In particular, 
it can be shown that from each simple zero of Q there emanate 3 Stokes lines along which 
Q{r)dr is imaginary. This means that one of the two solutions will grow exponentially 
whereas the other one will decay as we move away from the reference point t. On the other 
hand, we can also define anti-Stokes lines emanating from each zero of Q along which Q[r)dr 
is real and (f)f\ are oscillatory functions. By crossing anti-Stokes lines, the exponential 
behaviors of the two solutions are switched, while extending the solution across a Stokes 
line, the linear combination defining the solution is changed in a well-defined way: the 
coefficient of the dominant term stays the same, but that of the sub-dominant term picks 
up a contribution proportional to the coefficient of the dominant term. Thus we see that 
the appropriate contour that can be used will trace anti-Stokes lines and cross Stokes lines. 
Figure [7] shows such a contour for Schwarzschild geometry assuming (large) purely imaginary 
frequency. 

One can see from the potential ( |T3i) that for \'^uj\ oo, all zeroes of Q approach the 
origin and the potential may be approximated by its expansion around r = 0. In this 
case, the ^-dependent terms are negligible. They enter in subleading contributions Q, 
Near r = 0, keeping only the most singular term, the potential in (|T3|) . which describes 
integer-spin perturbations, becomes: 

^ 2{l+n) 

Hu.~^(2 + n)V-l)^ (45) 
and the resulting WKB equation is: 

ip"{r) + Ro{r)^{r) ~ (46) 

where 



Roir) = ( ^ j 
However, one can show 



- ((/ -!)(- + 2)^ + (n + l)(n + 3)) rf 



(47) 



43( 1 that the solutions so obtained by identifying Ro{r) with 
Q^(r), do not have the right behavior near the origin when compared with the exact solution 
near there. This can be rectified by using a slightly different R{r), which is permissible. The 
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modified form turns out to be: 



R{r) = i?o(r) - — ^ Q\r) 
For u = —i\uj\, the zeros of are: 



2(n+l) 



2 (n + 2)V 2{n+l)' 

Co" . . r, 



4^2(n+2) 



2|cj| ' ^ 



= 0,1,- ■■2(2 + n) - 1 



(48) 



(49) 



These zeros will serve as reference points for the phase integrals, because in doing so, the 
Stokes constants are either +i for traveling anti-clockwise (or —i for clockwise), and so, we 
will need to switch between reference points as we go around the contour: 



(50) 



with 



Qdr 



_ (n + 2)V 2(n+l)' 
4l^2^2(n+2) h 



1/2 



(51) 



where ?' and k label two consecutive zeros. Using the substitution y = r-^^^— ittt, followed 
hy y = coshx, we get: 



rk 



Qdr = ^ 



so that: 





1 ' 


1/2 


/ 












= ^^30 = 


~ u = — 






Y 



T 



(52) 



(53) 



It is also possible to consider a massless Dirac fermion propagating in the bulk, which can be 



relevant in split 



bulk fermions 



'ermion theories where rapid proton decay can suppressed by introducing 



48 



49( 1 that are superpartners to localized scalars and gauge bosons. A master 
equation potential based on a higher-dimensional Dirac equation potential was derived in 



Vv 



bulk-Dirac 



/ 



dr 



D-2 



+ f 



+ 



D-2 



(54) 



By expanding near r = 0, we find this potential to be proportional to 7'-(''+3n)/2^ which, for 
the present purpose, can be treated as zero when compared to the other bulk perturbation 
given by (1451) . At this approximation, therefore, we may use p = 1 for a bulk Dirac fermion. 
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Turning our attention to the actual calculation of grey-body factors, we see that since lo 
in the present case is not real, the conservation of flux will take on the following generalized 
form: 

T(cu)T(-cu) + R{uj)R{-uj) = 1 (55) 

and therefore, we consider two separate solutions representing the boundary condition at 
physical infinity; one with frequency u and the other with —uj. In the first case, the reflected 
wave at infinity is represented by (f)f^\r), while the ingoing wave by (f)2^\r). For — cj, the 
two solutions switch roles, but the calculation is effectively the same, and we shall combine 
them into one by writing the wavefunction in general as: 

= Xe^Vf'^ + l"e-*V^^ (56) 

where 6' is a phase angle. In the case of +uj, X = i?+ and Y = 1, whereas for —uj, X = 1 
and Y = R^. For ^{uM) > 0, the outgoing wave boundary condition at spatial infinity can 
be analytically continued to the anti-Stokes line labeled a in the figure. To find the the T's 
and i?'s we need to travel a path on anti-Stokes lines starting and ending at point a. The 
solution ipa at point a for ±uj is described by the wavefunction given in fl56l) with j = 1. 
Since we don't cross any Stokes lines in extending the solution to ti, the above expression 
remains valid there, but in going to point b we do cross a Stokes line, and to account for 
the Stokes phenomenon we need to add to ipa a contribution proportional to 02 with the 
constant of proportionality equal to the Stokes constant multiplied by the coefficient of the 
dominant function on the Stokes line, 0i in this case, so that: 

iJb = i^a- iXe'yt^ = Xe'yf'^ + [Ye-'' - iXe'') (57) 

To extend the solution to point c, we now need to change the lower limit of integration to 
have to as reference: 

^|Jb = Xe'^-'^^^f + {Ye-''^^"' - 1X6'^+'") 0^*°^ (58) 

In going to point c there will be no exchange in dominance, since we cross no anti-Stokes 
lines: 

V-c = ^6 - iXe''^-'''(pf^ = Xe'^-'''(j)'{'^ + {Ye"'^+"' - tXe'^^'"" - iXe'^-'") 0^*°^ (59) 
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We now go to point c', and in doing so we cross no Stokes lines, so that the combination 
f l59|) remains vahd, but now the integral is evaluated around a contour that loops around 
the pole r = Th- Replacing that contour by the one to the left of the pole, we have: 

where F is the integral encircling r = clockwise: 

Tbuik = j)Qdr = -2m Res Q = -ivr^l + (J^^ ~ -^i^^P (61) 

where /3 = 1/Th. To connect the solution to point d, we need to change the lower integration 
limit to t^: 

= Xe''e^^(t)^l'^ + {Ye''' - iXe''"^'" - iXe'') ^f^e"*^ (62) 

and now to connect to point c?, we need to cross the anti-Stokes line to get inside the loop, 
and this means that (pi and 02 will interchange dominance, further, the Stokes constant is 
now i, since we have reversed direction: 

= ^^, + i {Ye-'' - iXe''-^'^ - iXe'') ^^^e"*^ 

= [Xe^'e^^ + {lYe-'' + Xe^^'^- + Xe'') e'^^ (63) 
+ {Ye-'' - iXe"-^'^ - iXe") ^^^e"*^ 

Then switching back to point t^: 

ijj^ = IXe^'e'^e-'" + {iYe-^'e-'" + Xe''-^'" + Xe^'e-'") e"^^] 0f°^ 

(64) 

+ {Ye-^'e'" - iXe''-'" - iXe^'e'") 0^*°^e-^^ 
Going back to point c and crossing a Stokes line, we encircle the pole at the origin: 

^_ = [(2e-- + e'" + e-^'") e'^'-^^X + e*(-+r+e)x + i {e-'" + e'") e-^^^+^^F] 0?°^ 

+ e-*^ [-le'^-'^+'^X - le'^'^'^ + e'^""-'^] ^^^^ 

where the bar is to indicate that we have encircled the pole at the origin. Connecting to 
point h: 

^-^ = [(2e-'" + e"" + e-^'") e'^^X + e'^-'+^^X + i (e"^" + e^") e-'^^+^e)^^ ^ie^jto) 

+ {% {e-'^ e-'"''-'^ + e*^) e'^^'-^'^X - e-*{-+r+e)y] ^^^^ 
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(67) 



switching to reference point ti. 

ir, = [(2 + e^- + e-^'"") e'^^X + e'^X + t {l + e''") e-'^^+''^Y] e^Vf'^ 
+ (e-*r ^ ^-2^u-^^ ^ ^^^^^ e*(e-2.)^ _ g-.(2.+r+e)y] ^{ti) 

and finally returning to point a: 

^-^ = e-'^ \e'' (2(1 + cos(2z/)) + e'*^) X + 2ie'^''~'^ cos{u)Y] ^f^^ 

(68) 

+ e-'^ [Zie'^^-"^ cos z/ (l + e'^^ + 2 cos(2z/)) X - e-'^(l + 2 cos(2z/))r] 4*^^ 

To find R{uj) we need to properly impose the boundary conditions on X and Y, where Y will 
be set to unity, then we identify the coefficient of 0i in ipa as the same coefficient appearing 
in ipa except that because of the trip around the contour we have gained a phase: 

e-'^ [e'' (2(1 + cos(2z/)) + e'*^) R + 2te'^'-'^ cos(z/)] = e'^Xe''^ (69) 

Solving for R gives us the desired result: 

= -e^(-2.) ^^""^^ (70) 

^ ^ + 1 + 2 cos(2i^) ^ ^ 

To find the corresponding T, we notice that the boundary condition at the horizon, with (j)2 
dominating can be connected to point d, which was reached by moving along an anti-Stokes 
line along which 02 is dominant when u — > —ioo. Therefore, we can make the following 
identification: 

= ... + Te-*«0^)e-'r (71) 
comparing with the expression obtained in flM|) . we can solve for T{uj): 

T{uj) = -— (72) 

^ ' 1 + 2 cos(2z/) +ef^'^ ^ ^ 

For the case of —u, we just fiip the roles of X and Y, so that Y = R{—uj) is now obtained 
as follows: 

^- = . . . + i?(-cu)e-'V2*'^e'^ (73) 
with X = 1. By comparing with fl^S]) . we find: 

i?(-cj) = 2ie'(2e-.)pQg^ (74) 

Likewise, to find T{—uj), we make the identification 

= Te-^V^^e-^^ + ■ ■ ■ (75) 
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and by comparing with the expression in we find that: 



n-u) = 1 

The grey-body factor that results from this reads: 

7(0;) = T{uj)T{—uj) 



(76) 



1 + 2 cos(2z/) + e/^^ 



(77) 



This result is consistent with the asymptotic quasinormal frequencies found by Cho 4^ in 



4D, and it also implies Neitzke's results 



50| for scalar perturbation around a Schwarzschild 



hole for dimensions D > A, where the only modification due to extra dimensions appears 
through the Hawking temperature (cf. eq. ([2])). Further, in the 6D model that includes 
tension, one finds that an additional modification appears (cf. eq. (TT^). 

Moving on to brane-localized emissions, we need to use eq. (132|) as the radial master 
equation. This would describe scalar emissions {s = 0), massless Dirac fermions {s 



2n 



and gauge photons (s = 1). To apply our method, we need to put that equation in the 
Schrodinger form, in which case the potential looks like this 381]: 



VAr) 



f 



As + 9ns{l -f) + ^{in + 1)(1 - /) - 2f}' + sf 



lUJS 



(78) 



where 



and 



+ — {(n + l)(l-/)-2/} 
r 



qns = (2s + n + l)(ns + s + 1) - {s/2){n + l){n + 4) 



Ajs = 3{3 + l)-s{s + l) 

To also talk about the graviton localized to the brane, we have to use another potential [e.g 
3): 

VG{r) = ^ + l)-{{n + lf + 2} (1 - /)] (79) 

However, for our purpose, what we care about is the behavior of these potentials near r = 0. 
Keeping only the most singular term in both cases, we can combine the two potentials into 
one expression: 



Vbrane ^ (^r. 



r 



2(l+n) 

_h 

^2{2+n) 



(80) 
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with 

I \s\n-lf-\sn{n~l)~{n + l), fors = 0,i,l 
[ (n + 1)2 + 2, fors = 2 

To find the greybody factors in this case, we duplicate our steps for bulk emissions, and we 
find that the only difference occurs in the z/, which we can readily find by comparing the 
expressions flHUl) and fH5]) : 

t'brane = -TTW — -^^^ + ^ (82) 



In this case, the value of the contour integral Fbrane is found to be: 



= * = 1 - ' « _ 1 „^ (83) 

which is asymptotically the same as Fbuik- 

These results show that the grey-body factors in all cases are governed by (177|) with 
differences coming through the value of u. Tables H] and [III summarize this for bulk and 
brane-localized emissions, respectively. 

TABLE I: Expressions for v for different bulk perturbations 



Perturbation i^buik 
scalar &: gravi-tensor 
gravi-vector — tt 

EM vector Irrvr 

n+2 

EM scalar -^tt 

n+2 



gravi-scalar ~ -| [2 + 0.674(n + 4)"° 

Dirac fermion — ^ 



-0.54451 



Since the asymptotic quasinormal frequencies are the poles of the grey-body factors, we 
see that the standard formula still holds in form. 



27rz (n + i) +ln(l + 2cos2z/) (84) 
where n is the mode number, however, u now depends on the number of extra dimensions 



(as already found in 



5l|). 
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TABLE II: Expressions for v for different brane-localized perturbations 



Spin fbr 



" ~ 2(n+2) 

2 23/2 n+2 



^n(n-2)+2 
~^ 2(n+2) 



/ (n+l)2+2 1 
^ (n+2)2 + 4 



As a special case, the addition of tension in an n = 2 spacetime, will only modify Th, 
where it can be seen from equations (fT2l) and ([2]) that the Hawking temperature decreases 
with brane tension, 

Th ^ h^l^'T^ (85) 

Since this occurs in an exponent (eq. fl77j) ). even a moderate value of h will lead to 7(0;) — > 1 
faster with increasing frequency. 



V. CONCLUSION 



With the coming online of the LHC, the possibility of experimentally testing many of 
the theories that go beyond the Standard Model is a realistic one. Detecting black holes in 
high-energy collisions, if observed, would, for the first time, offer the opportunity to probe 
the realm of quantum gravity helping us resolve some of the most important puzzles of 20th 
century physics. This means that we need to be well-prepared with realistic predictions to 
interpret the possible observations. In the study of the creation and evaporation of mini 
black- holes, a main prediction of brane- world models, brane tension has often been neglected. 
While predictions may not be radically modified by including it, the fact that we are trying 
to observe something we have never observed before makes it imperative to strive to be as 
realistic and accurate as possible. Building on earlier work, we made an attempt to address 
this issue. 

We explored some of the implications of the codimension-2 model on grey-body factors 
for various types of bulk and brane emissions from a black-hole residing on a brane of finite 
tension. We saw that the results presented follow primarily from two modifications: (a) the 
enhancement of the horizon radius, which enhanced the emission spectra with increasing 
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tension, and (b) the increase of the angular eigenvalues with tension, which resulted in 
amplified absorption probability on the brane. We calculated analytically grey-body factors 
in the low frequency regime. We obtained expressions that had a non-trivial dependence 
on the brane tension with observable implications which we discussed. In particular, we 
saw that the observable emission rates got enhanced for both bulk and brane emission, yet 
we obtained a higher rate on the brane. It should be stressed that in the case of rotating 
black-holes our results are modified due to super-radiance. 

We also obtained analytic results in the case of large imaginary frequency. We derived 
expressions for grey-body factor in the bulk and on the brane for a Schwarzschild black hole 
in an arbitrary number of dimensions. The special case of two extra dimensions with tension 
was seen to follow from a simple modification of the Hawking temperature. 

We compared our analytic expressions in the two asymptotic regimes of low and high fre- 
quencies, respectively, with exact numerical solutions and discussed the range of intermediate 
frequencies where our analytic expressions are not accurate. The size of this intermediate 
range varies depending on the different parameters of the setup. Interestingly, this range is 
shifted toward high frequencies with increasing tension indicating a wider range of validity 
of our analytic low frequency expressions. 

In conclusion, our results show that brane tension has a pronounced effect on predictions. 
Although the model we discussed is not the most general one, it offers the opportunity to 
explore both qualitatively and quantitatively the effects of adding tension. It would be 
interesting to go beyond the model considered here generalizing it to other geometries and 
an arbitrary number of extra dimensions. We hope to report on progress in this direction 
soon. 
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FIG. 1: Low frequency absorption probability for bulk emissions as a function of b. Parameter 
p G {O, ^, |, 2, 2.25} from eq. increases from bottom to top. Here, £3 = 2, m = 1, = 0.2, 
and fi = 1. 
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FIG. 2: High frequency absorption probability for bulk emissions as a function of b. Parameter 
p G {O, ^, |, 2, 2.25} from eq. increases from bottom to top. Here, £3 = 2, m = 1, = 2.2, 
and = 1. 
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FIG. 3: Absorption probability for bulk scalars {p = 0) as a function of w for various values of 
brane tension, where h G {0.4, 0.6, 0.8, 1} increasing from bottom to top, with = 2, m = 1, /x = 1. 
The points on the graph come from exact numerical calculations. 
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FIG. 4: Bulk absorption probability for scalars for different values of m: for ^3 = 5, and m < £3, 
m increasing top to bottom, where b = 0.3, = 1. 



28 



I^branep 



10"^ 




FIG. 5: Absorption probability for brane-localized emissions as a function of b. The curves from 
top to bottom correspond to s = ^, s = and s = 1, respectively. Here, i = 2, lu = 0.02, and 
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FIG. 6: Absorption probability for brane- localized scalars (s = 0) as a function of uj for various 
values of brane tension, where b £ {0.4, 0.6, 0.8, 1} increasing from top to bottom, with i = 2,m = 
l,fi = 1. The points on the graph come from exact numerical calculations. 
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FIG. 7: Stokes lines structure for the Schwarzschild black hole in 6D in the complex r-plane for 
large imaginary frequency. Empty circles are the zeros of Q and the filled circle is the horizon 
radius. 
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FIG. 8: Energy emission spectrum for scalars in the bulk for various values of brane tension, with 
h £ {0.4, 0.6, 0.8, 1} increasing from top to bottom. 
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FIG. 9: Energy emission spectrum for scalars on the brane for various values of brane tension, with 
b £ {0.4, 0.6, 0.8, 1} increasing from top to bottom. 
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FIG. 10: Relative bulk-to-brane emissivity for scalars for different values of brane tension, with 
b G {0.4, 0.6, 0.8, 1} increasing from top to bottom. 
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FIG. 11: Energy emission spectrum for gauge vectors in the bulk for various values of brane tension, 
with h G {0.4, 0.6, 0.8, 1} increasing from top to bottom. 
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FIG. 12: Energy emission spectrum for gauge vectors on the brane for various values of brane 
tension, with b G {0.4, 0.6, 0.8, 1} increasing from top to bottom. 
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FIG. 13: Relative bulk-to-brane emissivity for gauge vectors for different values of brane tension, 
with b S {0.4, 0.6, 0.8, 1} increasing from top to bottom. 
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